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Abstract 

The B c meson is the bound state of be (or bc) whose recent detection 
is the first step toward completion of the spectroscopy of heavy quark 
mesonic states. The b-c states have properties that conveniently fill 
the gap between the J/ij) and the T states. Thus it is probable that at 
RHIC the B c mesons will serve as a probe of deconfined matter. We 
find that significant differences arise for B c formation in deconfined 
and confined matter. Our initial calculations suggest that: 

(a) The rates of normal hadronic production mechanisms at RHIC 
energies are not sufficient to produce a detectable number of B c mesons. 

(b) If a region of deconfined quarks and gluons is formed, the 
production (and survival) rate can be enhanced by several orders of 
magnitude. 

(c) The observation of B c mesons at RHIC would signal a source 
of deconfined charmed quarks, and the rate of B c production will be 
a measure of the initial density and temperature of that source. 



1 Introduction 

This work||U] investigates the possibility that the production of B c mesons at 
RHIC may serve as a signal for the presence (or absence) of a deconfined state 
of matter 0. The study of the b-c sector has the advantage of a long history 
of potential model analysis in the bb and cc sectors. These studies have 
provided robust predictions for the mass and lifetime of the B c states 0, and 
the recent measurements by CDF|4| are consistent with those calculations. 

First let us estimate at RHIC the production rate of different heavy quarks 
and mesons, which one would expect if it results just from a superposition 
of the initial nucleon-nucleon collisions. For heavy quark production, pQCD 
calculations for p-p interactions fit present accelerator data and bracket the 
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RHIC energy range. Hard Probes Collaboration || estimates indicate about 
10 cc pairs and 0.05 bb pairs per central collision at RHIC. J/ijj and T pro- 
duction involves the use of some model, such as the Hard Probes color singlet 
fits||, which would predict bound state fractions of order somewhat less than 
the one percent level. 

A similar analysis for B c production reaches significantly different results^ 
Since the b and c must be produced in the same nucleon-nucleon interaction, 
parton subprocesses of order a A s are the leading order contributions. This 
leads to a substantial reduction of the bound state fraction 

_ B c + Bl 

relative to the few percent levels for the corresponding T state fractions. At 
RHIC energies, typical values are Rb = 3 — 10 x 10~ 5 , with the uncertainty 
from the scale choice in the pQCD calculations 0. 

To convert these numbers into B c production predictions for RHIC, we 
have looked at two scenarios for the luminosity, a) The "first year" case as- 
sumes a luminosity of 20 inverse microbarns with no trigger, b) The "design" 
luminosity assumes 65 Hz event rate with a 10% centrality trigger in Phenix, 
and uses 10 7 sec/year. The predictions we obtain are listed in Table 1. In- 
cluded in the estimates are both the weak branching fraction of the B c plus 
the dimuon decay fraction for J/ip. Similar numbers are shown for the J/i/j 
and T production and detection via and also the underlying heavy 

quark production which may be useful to make contact with other estimates. 
One sees easily that in this scenario there is no hope of seeing 5 c 's 
at RHIC. 



Table 1: RHIC yields for heavy quark systems. 



Observable events 


First Year 


Design Luminosity 


cc-pairs 


2.810 s 


6.5 10 9 


66-pairs 


1.2 10 6 


3.2 10 7 




1.6 10 5 


3.9 10 6 


T(ls) — > fi + ^r 


140 


3800 


D 2.5% T/ii 6% — i— — t 




(No Deconfined Phase) 


0.05-0.18 


1.5-4.9 


(QGP+cc in Chemical Equil.) 


18 


490 


(Only initial cc at T = 500 MeV) 


130 


3530 


(Only initial cc at T = 400 MeV) 


235 


6420 


(Only initial cc at T D = 300 MeV) 


475 


12900 



1 In the following we include in the term B c also the vector IS state B*, since its mass 
splitting should only allow an electromagnetic decay into the pseudoscalar ground state 
and thus both will contribute identically in the experimental signatures. 
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2 Deconfinement Scenario 

Now the principal reason for our interest - could deconfinement change the 
B c production rate at RHIC? We have investigated the following scenario: 
In those events in which a bb pair are produced, one could avoid the small 
B c formation fraction if the 6-quarks are allowed to form bound states by 
combining with c-quarks from among the 10 cc pairs already produced by 
independent nucleon-nucleon collisions in the same event. This can occur if 
and only if there is a region of deconfinement which allows a spatial overlap 
of the b and c quarks. In addition, one would expect some cc production in 
the deconfined phase during its lifetime, as a result of the approach toward 
chemical equilibration. The large binding energy of B c (840 Mev) would 
favor their early "freezing out" and they will tend to survive as the temper- 
ature drops to the phase transition value. The same effect for the B mesons 
and indeed for the B s will not be so competitive, since these states are not 
bound at the initial high temperatures (or equivalently they are ionized at a 
relatively high rate by thermal gluons). 

To do a quantitative estimate of these effects, we calculate the dissociation 
rate of bound states due to collisions with gluons, utilizing a quarkonium 
break-up cross section based on the operator product expansion ||: 



where k is the gluon momentum, e Q the binding energy, and /i the reduced 
mass of the quarkonium system. This form assumes the quarkonium system 
has a spatial size small compared with the inverse of Aqcd, and its bound 
state spectrum is close to that in a nonrelativistic Coulomb potential. The 
magnitude of the cross section is controlled just by the geometric factor -p-, 
and its rate of increase in the region just above threshold is due to phase 
space and the p-wave color dipole interaction. 

For the breakup rate A# of B c states in deconfined matter, we calculate 
the thermal average: 



where v g = 1 and all modes of massless color octet gluons have been included. 
Numerical results for these rates are shown in Fig. [|. For comparison, 
breakup rates are also shown for the J/ip and T (and even the B s , but 
the the approximations made for this cross section probably have a very 
marginal validity in view of such a large state). One sees that in the range 
of temperatures expected at RHIC, these breakup rates for B c lead to time 
scales of order 1 — 10 fm. 




(1) 




(2) 
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Figure 1: Thermal QGP quarkonium dissociation rates as functions of tem- 
perature. 

For an estimate of the corresponding cross section for the formation re- 
action a Fib + c — > B c + g) we utilized detailed balance relations. This leads 
to a finite value of op at threshold, since it is an exothermic reaction. In the 
approximation that the massive 6-quarks are stationary, which is expected 
to be a reasonable approximation due to their energy loss in the hot plasma 
f|, the formation rate is then calculated for a thermal distribution of charm 
quarks: 

\ F = (v c n c a F ) = — tt P 2 dp e~~ a F {p) (3) 



where E p = yp 2 + m? c . These formation rates are shown in Fig. They 
have been calculated for three different values of charm quark mass. It is 
apparent that the results are quite sensitive to this choice, due to the strong 
dependence of total charm quark population. The same values of m c have 
very little effect on the breakup rates, since they only change the overall scale 
in the geometric factor of the breakup cross section. 
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Figure 2: Thermal QGP B c formation and dissociation rates as functions of 
temperature. 



Also shown in Fig. are the ratios Xs/Xp, which in our normalization is 
related to the bound state fraction in the equilibrium limit:[| 



3 \ F 



B c + B* 2A B (A s 

f " = -ffl— TTjg- (4) 



Note that this ratio approaches its upper limit of 2 when the formation rate 
dominates over the breakup rate. This corresponds to the situation in which 
every b-quark produced in the initial collisions emerges as a B c bound state. 

We choose a transition temperature Tf = 160 MeV at which to evaluate 
the final bound state populations. Here the equilibrium bound state fraction 
Rb drops to as low as several percent, but it is at least a factor of 100 above 
what one may expect in the no-deconfinement scenario. We have chosen to 
use the equilibrium ratios although at this final temperature the rates are not 
sufficient for them to be approached. This provides an even more conservative 
estimate for the final bound state populations. The corresponding entries in 
the Table for numbers of B c mesons (labeled QGP + cc in Chemical Equil.) 
uses this conservative lower limit estimate. 

Implicitly, this analysis uses the full chemical equilibrium density for c- 
quarks. To get a more realistic limit we repeated the calculation, using 

2 This bound state fraction is reached if the system has enough time in its dynamical 
evolution to relax to the steady-state solution at each temperature. We have verified that 
this is roughly the case down to T = 300 MeV, at which point the B c abundance begins 
to freeze out. 
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only the initially-produced c-quarks in the formation rate. From the initial 
population of 10 cc-pairs produced via nucleon-nucleon collisions in a central 
Au-Au collision at RHIC, and an initial volume V Q = tt(Rau) 2 ^o with r G = 
1.0 fm, one concludes that only for initial temperatures T Q < 300 MeV is the 
initial charm quark density comparable to that for full chemical equilibrium. 
For initial temperature T = 500 MeV, for example, the chemical equilibrium 
charm quark density would be about a factor of 40 higher than that actually 
provided by the initially-produced charm quarks. As temperature decreases 
below T , the isentropic expansion VT 3 = Const, leads to a decrease in 
the c-quark density proportional to T 3 , rather than the e~ mc l T of chemical 
equilibrium. We have verified that the rates of both charm annihilation and 
production in a deconfined state for T < 300 MeV then lead to charm quark 
occupancies which exceed those for chemical equilibrium as one approaches 
the transition point [10]. Fig. |3] displays a comparison of chemical equilibrium 



charm quark densities and those resulting from a constant number of initially- 
produced charm quarks with isentropic expansion. 



10 u 



E 10" 



10" 



N c = 1 0, m c = 1 .5 GeV, x = 1 fm, VT = const. 



T = 500 MeV 
T = 400 MeV 
T = 300 MeV 

Chemical Equilibrium Density 




Figure 3: c-quark density from initial production at RHIC. 



These more realistic charm quark densities are used to recalculate the 
formation rates, and the resulting ratios Xf/^b are shown in Fig. |] for several 
values of initial temperature T Q . The last few rows in the Table show the 
corresponding B c numbers at RHIC in this scenario, where we have used the 
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equilibrium bound state fractions again at a final temperature Tf = 160 MeV. 
They depend quite strongly on the initial temperature, which determines the 
final charm density through the assumed isentropic expansion. 




Figure 4: Ratio of formation to break-up rates of B c as function of temper- 
ature for fixed charm quark abundance. 



We are in the process of refining these preliminary results JTOJ. Initial 
numerical solutions of the kinetic equations using time- dependent formation 
and breakup rates indicate the final bound state populations saturate at val- 
ues appropriate to those for equilibrium temperatures somewhat above the 
transition values. This would be expected, since the rates at low tempera- 
tures are not sufficient to reach the equilibrium solutions before the volume 
expansion reduces the temperature to even lower values. Also, production 
and annihilation of additional charm quark pairs is most effective at higher 
temperatures, which enhances the effective formation rates. Both of these 
effects will enhance the bound state production fractions for the higher initial 
temperatures, and reduce it somewhat for lower initial temperatures. How- 
ever, it appears that the sensitivity to the parameters of the deconfined state 
will remain, making the B c signal a sensitive probe of QGP. 

While numerical considerations presented here will see a considerable re- 



finement in the near future fllP] , the firm conclusion we are able to make 
today is that should QGP be formed at RHIC there would be a very sig- 
nificant enhancement of the formation of B r mesons which can be observed. 
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The primary mechanism responsible for this enhancement is the interaction 
of initially-produced bottom and charmed quarks, which will not operate in a 
confining phase. The observation of any I? c 's at RHIC is thus both a "smok- 
ing gun" signal of deconfinement and a probe of the initial temperature of 
the system and the initial density of deconfined charm. 
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